Adiabatic dynamics of periodic waves in Bose-Einstein condensate with time 

dependent atomic scattering length 
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Evolution of periodic matter waves in one-dimensional Bose-Einstein condensates with time de- 
pendent scattering length is described. It is shown that variation of the effective nonlinearity is a 
powerful tool for controlled generation of bright and dark solitons starting with periodic waves. 
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Observation of Bose-Einstein condensate (BEG) in 
gases of weakly interacting alkali metals have stimulated 
intensive studies of the nonlinear matter waves. A new 
area of physics — nonlinear matter waves and nonlinear 
atomic optics — was originated. Generation and dynam- 
ics of solitary wave pulses in BEG's is one of the most 
important related problems. Experimental observations 
of dark and bright 0,0 solitons have recently been 
reported. Theoretically, several methods of creating soli- 
tary waves have been proposed. First of all, it is a mod- 
ulational instability , which is a universal phenomenon 
of the nonlinear physics (especially intensively explored 
in nonlinear optics, see This method, however, can- 
not predict exactly the parameters of generated solitons. 
Another method, which is controllable in the above sense, 
is the so-called phase engineering j2i], which consists of 
imposing an initial phase on a BEG and is appropriate for 
generating dark solitons. The phase imprinting, however, 
affects the whole background condensate which acquires 
nonzero initial velocity and starts to oscillate in a trap 
potential. The problem becomes even more complicated 
when one is interested in generating trains (or lattices) 
of solitons in BEG's. 

In the present Letter we show that a powerful tool for 
generating and managing matter soliton trains can be 
provided by variation of the effective nonlinearity, which 
in practical terms can be achieved by means of variation 
of the s-wave scattering length due to the Feshbach 
resonance 0: 



asit) = a{l + A/iBo-B{t)). 



(1) 



Here a is the asymptotic value of the scattering length 
far from resonance, B{t) is the time-dependent external 
magnetic field, A is the width of the resonance and Bq is 
the resonant value of the magnetic field. Feshbach reso- 
nances have been observed in Na at 853 and 907 G Q Lin 
^Li at 725 G H, and in ^SRb at 164 G with A = IIG 
Also, rapid variation in time of as has been recently used 
for generation of bright solitons in BEG 0, 0| . Here we 
want to indicate that in quasi-one-dimensional geometry 



an initially weak modulation of the condensate can be 
amplified by means of proper variation of the scattering 
length. As a result, the condensate evolves into either 
a sequence of bright solitons for < 0, or "domains" 
separated by dark solitons for as > 0. In the case of 
bright solitons the attractive forces between atoms ex- 
actly compensate the wave-packet dispersion in the lon- 
gitudinal direction, so that the confining trap potential 
in this direction becomes unnecessary. Then the motion 
of bright solitons in the longitudinal direction can be con- 
trolled by means of application of external forces. Actu- 
ally, oscillations of bright solitons in the trap observed in 

Q give a simple example of such controllable motion. 
Thus, quasi-one-dimcnsional bright BEG solitons behave 
as separate entities and their investigation seems to be 
a quite promising field of research. Dark solitons may 
be considered as moving "domain walls" which separate 
regions of a condensate with different values of the order 
parameter. Investigation of dark solitons is also useful 
for understanding the properties of BEG. In general, the 
problem of the controllable soliton generation is impor- 
tant for a number of BEG applications, like atomic inter- 
ferometry llOl. and different kinds of the quantum phase 
transitions 11], as well as in the context of the nonlinear 
physics, including nonlinear optics and hydrodynamics. 

Our approach is based on the well established con- 
cept that the BEG dynamics at low enough temperature 
is well described by the three-dimensional (3D) Gross- 
Pitaevskii (GP) equation. In some physically important 
cases it admits a self-consistent reduction to the ID non- 
linear Schrodinger (NLS) equation 
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In particular this is the case of a cigar-shaped BEG of 
low density when e = ^^^^^ <^ 1 and ^ 
is a total number of atoms, a± — { — 
( — - — )^/^ are linear oscillator lengths in the transverse 
and in cigar-axis directions, respectively (in the small am- 
plitude limit they are of order of effective sizes of the con- 
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densate) , wj^ and ujq being respective harmonic oscillator 
frequencies, < 1 is a positive parameter , and nia is the 
atomic mass. In Q time t and coordinate x are mea- 
sured in units 2/{e^Lu_i) and aj_/€ respectively. The order 
parameter in the leading order is searched in the form 

where a = sign(a5), and g(t) = as{t)/as{0). It will be 
assumed that ae{t) does not change its sign and thus g{t) 
is a positive- valued function. We notice that smallness of 
the density rules out a possibility of collapse phenomenon 
(if a < 0). 

We start with analytical estimates supposing that the 
initial wave function u(x, 0) is modulated along the cigar 
axis with the wavelength L of modulation much less than 
the longitudinal dimension of the condensate, i.e. of 
the I: L <^ I. Therefore, at this stage we neglect the 
smooth trap potential and impose cyclic boundary con- 
ditions. Then the initial wave function can be approxi- 
mated well enough by exact periodic solutions of Eq. (0) 
a.t v — 0. For example, if at < = we take into ac- 
count only one space harmonic of the initial wave func- 
tion, u{x,0) — uq -f Ml cos(x/L), then this distribution 
can be approximated by well-known elliptic function so- 
lutions of Eq. Q with a small parameter m (see below) . 
We are interested in evolution of such solutions due to 
slow change of g{t) with time. At the same time we sup- 
pose that the total time of adiabatically slow change of 
g{t) is much less than the period ^ 2-k/v oi soliton oscil- 
lations in the trap potential, so that we can neglect the 
influence of the trap potential on the motion of solitons 
during the formation of soliton trains and put = in 
Eq. ||2Jl. This means that we shall consider analytically 
relatively small segments of the modulated wave much 
greater than the wavelength L and much smaller than 
the size I of the whole condensate in the trap. To solve 
the problem of the condensate evolution, we note that 
substitution 

u{x,t)^v{x,t)/^/^ (3) 

transforms Eq. ||2Jl with v = {] into 

ivt + Vxx — 2a\v'\^v — iev (4) 

with e{t) = g'{t)/2g{t). Thus, for slowly varying g{t) the 
right hand side of Eq. Q can be considered as a small 
perturbation: |e(i)| <C 1. As it follows from Eq. 
for the initial distribution one has ■i;(a;, 0) = u(a;,0). For 
our purposes it is enough to consider typical particular 
solutions of the unperturbed NLS e quat ion which are pa- 
rameterized by two parameters Ai_2 uM- Under influence 
of the perturbation, these parameters in the adiabatic ap- 
proximation become slow functions of time t. Equations 
which govern their evolution can be derived by the fol- 
lowing simple method. 

First, the initial values of Ai_2, as well as the coefficients 
in Eq. are supposed to be independent of x, hence the 



wavelength L of the nonlinear wave evolving according to 
Eq. is constant, 

dL{\i{t),\2{t))/dt^Q. (5) 

Second, we find that the variable N {\i{t) , X2{t)) — 
jwpdx changes with time according to 

diV(Ai(i), \2{t))/dt = 2eN{\i{t), X^it)). (6) 

Then, if the expressions for L and N in terms of Ai^2 
are known, Eqs. |(3J| and reduce to two equations 
linear with respect to derivatives dXi 2/dt, which yield 
the desired system of differential equations for Ai_2- The 
form of this system depends, of course, on the choice 
of the parameters Ai^2- It is well known from the the- 
ory of modulations of nonlinear periodic waves that for 
completely integrable equations (as the NLS equation) 
the most convenient choice is provided by the so-called 
"finite-gap integration method" of obtaining periodic so- 
lutions. Therefore we shall use the parametrization of the 
periodic solutions of the NLS equation obtained by this 
method (see, e.g. 0|), and consider three most typical 
cases. 

Case 1: cn-wave in a BEC with a negative scattering 
length. In the case of a BEC with negative scattering 
length, a = —1, there are two simple two-parametric 
periodic solutions of unperturbed Eq. One of them 
has the form 

v{x,t) = 2Aie-4^(^i-^')*cn[2^A2 + A2x,m], (7) 

where the parameter of elliptic function is given by m = 
Xf/{Xi -I- A2). Then straightforward calculations give 

L^^m^' N = 8^Xl + XlEim)^AXlL (8) 
V Aj + A2 

where K(to) and E(m) are complete elliptic integrals of 
the first and the second kind, respectively. Substitution 
of these expressions into Eqs. © and © yields the sys- 
tem 

dXi _ ((A^ + A|)E(m) - A^K(TO))E(m)Ai 

"d7 ~ AfE2(TO) + A2(K(m)-E(m))2 ' 

dM _ (AiK(m) - (Xl + A^)E(TO))(K(m) - E(m))A2 

"d7 ~ A?E2(to) + A^(K(to) - E(to))2 

(9) 

where 

z^z{t)^2[ e{t')dt' ^\ng{t), z(0) = 0. (10) 

If dependence of Ai and A2 on z is found from ((SJ, then 
Eq. gives evolution of the periodic wave u{x,t) with 
slow change of the parameter z connected with time t by 
Eq. (|10|l . In particular, the density of particles in the 
condensate is given by 

|u|2 = Ae-'Xl{z) cn^[2^Jxl{z) + A^z) x, m], (11) 
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FIG. 1: Numerical solution of Eq.© with g{t) = e*''^. Ini- 

tial conditions are chosen in the form u{x^ 0)6'"" where 
u{x,0) is given by (a) Eq. Q with Ai(0) = 1, A2(0) = .2, 
T « 2 1/ = .02; (b) Eq. with Ai(0) = 1.5, A2(0) = .2, 

T 1.5 and u = .01; (c) Eq. (O with Ai(0) = 3, A2(0) = .3, 
r ~ 1 and u = .01. In the boxes we show density distributions 
at initial (thin lines) and final (thick lines) moments of time. 

where transformation to the time variable should be per- 
formed with the use of Eq. ifTUIl . 

In Fig. ^ we present an example of the evolution of 
the respective density distribution in the presence of a 
harmonic trap potential where the parabolic parameter 
as well as experimentally feasible parameters are used. 
The figure shows that in the case of a negative scattering 
length given by (QJ increase of the magnetic field B{t) 
within the region B{0) < B{t) < Bq results in compres- 
sion of the atomic density and formation of a lattice of 
matter solitons. 

Case 2: dn-wave in a BEC with a negative scattering 
length. Another simple solution of the NLS equation Q 
with (T = — 1 is given by 

v{x, t) = (Ai + A2)e-2»(^i+^')*dn[(Ai + Aa)^, m], (12) 

where m = 4AiA2/(Ai + ■ By analogy with (O we 
derive the following equations for Ai and A2: 

d\i _ Ai(Ai + A2)E(m) 

dz " (Ai - A2)K(m) (Ai A2)E(m) ' .-.o^ 
dA2 ^ A2(Ai -|-A2)E(m) ^^"^^ 

dz ~ (Ai-A2)K(m)-(Ai+A2)E(m)' 

where it is supposed that Ai > A2 and z is defined by 
Eq. H10() . Now the density of particles is given by 

|up = e-'{\i{z) + \2{z)fdn\\i{z) + \2{z))x,m]. (14) 



An example of the respective evolution in the presence 
of the potential is given in Fig. ^p. One again observes 
formation of a lattice of matter solitons starting with a 
weakly modulated periodic wave. 

Case 3: sn-wave in a BEC with a positive scattering 
length. In the case of (0)) with a — 1 there exists simple 
periodic solution 

v{x, t) = (Ai - A2)e2'(^?+^')*sn[(Ai + A2)a;, m], (15) 

where ni — and it is supposed that Ai > A2. 

Now we obtain 

d\i _ Ai(Ai + A2)(K(m) - E(m)) 

~ih ~ 2AiK(to) - (Ai + A2)E(m) ' 

dA2 ^ A2(Ai +A2)(K(m) -E(m)) ^^^^ 

dz ~ 2A2K(to) - (Ai + A2)E(m) ' 

with z defined again by Eq. 1)1 0|l . The density of particles 
in the condensate is given by 

\u\^ - e-^(Ai(z)-A2(z))W[(Ai(z)-HA2(z))a;,m]. (17) 

This case, but in the presence of external trap potential 
is illustrated in Fig. where by means of increase of 
the magnetic field a periodic wave is transformed into a 
lattice of dark solitons. 

In physical units, the cases depicted in Fig. ^ corre- 
spond to {&) N — 1.4 • 10^ Li atoms in a trap with 
aj^ « 7^m, and oq ~ 230 /um, ih) N = 2- 10"* Li atoms 
in a trap with a_L « 6 ^m, and oq ~ 416 /im, and (c) 
M = 10"* "^^Na atoms in a trap with a± = 3.4 /xm, and 
ao = 264 fim. In the last case, however one observes shifts 
of the soliton positions as well as decrease of a density of 
particles located about the potential minimum because of 
weak oscillations of the condensate in the trap potential 
(the expanding phase is depicted in the figure). Thus, 
although effective r corresponding to physical time 
to = 40 ms (used in all simulations) is not large charac- 
teristic amplitudes of solitons placed in the center of the 
trap potential (in the cases (a) and (b) match well the 
amplitude values following from the adiabatic approxi- 
mation developed for a homogeneous NLS and for large 
T, and in the case (c) one observes qualitative agreement. 
No instabilities of periodic waves are observed during pe- 
riods of soliton train formation. 

The developed analytical approach can be generalized 
to the NLS equation with linear damping, when the right 
hand side of Eq. ^ is equal to —iju, 7 being the damp- 
ing constant, the substitution Q yields again Eq. (0J) 
but now with modified value of e: e — > e — 7. Hence, the 
equations for Ai and A2 hold their form with z{t) defined 
as z{t) = In g{t) — 2jt, so that the only change in Eqs. 
(|ll|l . 1)141) . and 1)171) is multiplication of their right hand 
sides by exp(— 27^). Another generalization corresponds 
to moving soliton trains what may be useful for treatment 
of BEC in ring traps. The above method remains actu- 
ally the same, but the initial condition includes a phase 
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factor linearly depending on the space coordinate. Then 
the results presented in Fig. ^ can be viewed as plots of 
the respective current as functions of time. In practice, 
solitons can be put also in motion by imposing a proper 
external potential. One of applications of such moving 
soliton trains could be a "laser" of matter solitons. 

The consideration provided above implies that one 
starts the adiabatic deformation with an initially periodic 
solution. A question arises about possibility of creation 
of such a state experimentally. A natural approach to 
solving this problem would be the use of an optical trap 
01 . In such a trap it is possible to create a nonlinear pe- 
riodic distribution of a BEG 0, 0| . Then, switching off 
the laser beams, producing the trap, will result in a pe- 
riodic distribution of the condensate. However, it is not 
stable without the trap since it is not a solution of the 
respective GP equation. This difficulty can be overcome 
if simultaneously with switching off the optical trap one 
abruptly changes the scattering length (or alternatively 
provides change of the number of particles) in such a way 
that the existing distribution will satisfy To be more 
specific, let us consider an example of a BEG with a pos- 
itive scattering length in an optical trap given by [15|: 
V{x) = — 2Vosn(2Ka;, to), where Vq is the potential am- 
plitude, K = Xi + X2 and TO is the same as in (jl7|l . The 
equation describing BEG evolution now admits a solution 

u{x,t) = ^Vo + (Ai - A2)2 e2*(^i+^^)*sn[(Ai + \2)x, to]. 

This last function solves also Q with a = 1 and 
g — and thus by switching off the potential 

V{x) with simultaneous changing the scattering length 
by Afls = (Ar-A2)^ ^^"^ achieves the desired initial state. 
Notice that although experimentally a sn-potential is not 
easily reachable in a general case, for a large range of m it 
is approximated very well by only a few first Fourier har- 
monics. For example, for a situation described in Fig. ^ 
one has V{x) « -2Fo[1.47 sin(0.78x) + 0.15 sin(2.3x)] 
with the accuracy about 1%. 

In order to estimate characteristic scales of adiabatic 
deformations we introduce an "aspect ratio" defined as 
6 — \ AA\/L, where AA is a total variation of the ampli- 
tude of the periodic wave and L is its wavelength. Then 
the cases 6^1 and (5 <C 1 correspond to a solitonic lat- 
tice and to a modulated plane wave. Dependence of the 
scattering length on time can be simulated by g{t) = e*/*" 
(physical units). Taking as an example the solution de- 
picted in Fig. where (5 « 1 at t = 0, we find that al- 
ready at t « 30ms the aspect ratio becomes 6 « 100. The 
adiabaticity of the process means here that ^ ^ ^ ^ ^ 
(physical units). It is satisfied well enough for traps with 
w_L > 27r X 200 Hz and luq < 2'k x 5 Hz. We notice that 
the imposed condition rules out a possibility of resonant 
phenomena which might be related to coincidence of a 
linear oscillator period, originated by one of the trap di- 
mensions, on the one hand and the characteristic time of 



change of Us on the other hand. Stability of the above 
solutions has been studied numerically in where it 
has been found that soliton trains are stable in the case of 
positive scattering length and are also stable in the case 
of negative scattering length for special choice of the pa- 
rameters. In this context, the adiabatic variation of the 
scattering length, which results in the change of the wave 
parameters, can be used for stabilizing (or destabilizing) 
respective periodic solutions. 

To conclude, we have outlined the main idea of man- 
agement of periodic nonlinear waves in BEG's. The the- 
ory, although being developed for a homogeneous NLS 
equation, give an accurate estimate of the central part of 
a BEG placed in a magnetic trap, the later being studied 
numerically. The existence of trap can also be taken into 
account in the framework of more sophisticated theory 
recently developed in ^3 • Application of this method to 
the theory of a BEG will be presented elsewhere. 

Work of A.M.K. in Lisbon has been supported 
by the Senior NATO Fellowship. A.M.K. thanks 
also RFBR (Grant 01-01-00696) for partial sup- 
port. V.V.K. acknowledges support from the Eu- 
ropean grant, GOSYG n.o. HPRN-GT-2000-00158. 
Work of V.A.B. has been supported by the FGT 
fellowship SFRH/BPD/5632/2001. Gooperative work 
has been supported by the NATO-Linkage grant No. 
PST.GLG.978177. 



[1] S. Burger, et. al. Phys. Rev. Lett. 83, 5198 (1999); 
[2] J. Denschlag et al. Science 287, 97 (2000). 
[3] L. Khaykovich et al. Science, 296, 1290 (2002). 
[4] K.E. Strecker et al. Nature, 417, 150 (2002). 
[5] V.V. Konotop and M. Salerno, Phys. Rev. A 65, 021602 
(2002). 

[6] F.Kh. Abdullaev, S.A. Darmanyan and J. Garnier, In 

Progress m Optics, Edit. E. Wolf, 44, 303 (2002). 
[7] A.J. Moerdijk, B.J. Verhaar, and A. Axelsson, Phys. Rev. 

A 51, 4852 (1995). 
[8] S. Inouye, et al., Nature (London), 392, 151 (1998). 
[9] P. Courteille et al., Phys. Rev. Lett., 81, 69 (1998). 
[10] B.P. Anderson, K. Dholakia and E.M. Wright, eprint; 

cond-mat/0209252 
[11] R. Kanamoto, H. Saito, and M. Ueda, eprint: 

cond-mat/02 10229 
[12] Our choice of the parameters is dictated by simplicity 

of the obtained equations. Physical parameters can be 

expressed in terms of Ai (see lO, 1121 1 or 115l l'). 
[13] A.M. Kamchatnov, Nonlinear Periodic Waves and Their 

Modulations — An Introductory Course, (World Scientific, 

Singapore, 2000) 
[14] B.P. Anderson and M. Kasevich, Science 282, 1686 

(1998); S. Burger et al., Phys. Rev. Lett. 86, 4447 (2001); 

O. Morsch, et al., Phys. Rev. Lett. 87, 140402 (2001). 
[15] J.C. Bronski, et al., Phys. Rev. Lett., 86, 1402 (2001). 
[16] L. D. Carr J. N. Kutz, and W. P. Reinhardt, Phys. Rev. 

E 63, 066604 (2001). 
[17] A.M. Kamchatnov, eprint: |nlin.SI/0302027l 



